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Discriminating between ^ Vx and ^ ^ sterile in atmospheric v 
oscillations with the Super-Kamiokande detector. 
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Abstract. A strong body of evidence now exists for atmo- 
spheric V;j disappearance oscillations. Such disappearance 
could be explained by oscillations to either Vx or a "sterile" 
neutrino (Vy). Super-Kamiokande uses three different meth- 
ods to distinguish between these two scenarios. First, matter 
effects would suppress the ^ V.s oscillation amplitude at 
high energy. Second, oscillation to v , would reduce the over- 
all neutral-current neutrino interaction rate. Third, the smok- 
ing gun of ^ Vx oscillations would be the observation of 
1 appearance resulting from charged-current Vx interactions. 
The results of these three techniques are presented, which 
strongly favor ^ Vx oscillations over v^ ^ Vj. 



1 Introduction 

The Super-Kamiokande (Super-K) experiment has reported 
evidence for the oscillation of muon neutrinos produced in 
cosmic -ray induced showers in the atmosphere, via the ob- 
servation of V^ disappearance as a function of neutrino path- 
length and energy (Y. Fukuda et al, 1998, 1999). The ob- 
served Vf, signal shows no excess, so is not consistent with 
substantial v^j <-> Vg oscillation; therefore the missing v^ must 
be changing to an "invisible" flavor. Since a Vx charged- 
current (CC) interaction would produce a i lepton, Vx below 
the 3.4 GeV i production threshold do not interact via the CC 
channel. Vx above this threshold would produce T leptons, 
but 1 decay products would produce a multiplicity of parti- 
cles. Since the data set used to observe v^j oscillation selects 
only single-particle events, only the resulting disappearance 
of V^ is thus observed. 

An alternative scenario that can explain muon neutrino dis- 
appearance is oscillation with a sterile neutrino (Vj), which 
undergoes neither CC nor neutral current (NC) interactions. 
A fourth flavor of neutrino must be invoked to simultane- 
ously explain apparent neutrino oscillations in three widely 
separated regimes: atmospheric neutrinos; solar neutrinos 
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(Bahcall et al. , 1998); and flie LSND experiment (C. Athanas- 
sopoulos et al., 1998). However, studies of the Z" at LEP 
(Decamp et al, 1992) indicate that any such new light neu- 
trino flavor must not interact via the weak force, leading to 
speculation about a sterile flavor 

Super-K has used two methods to distinguish between the 
Vx and Vj hypotheses for explaining atmospheric v^ disap- 
pearance (Y. Fukuda et al., 2000). The first method is to ex- 
amine events likely to have been caused by NC interactions. 
While Vx readily undergo such interactions, Vj would not, re- 
sulting in a relative suppression of the NC signal. The second 
method exploits so-called "matter effects", which reduce the 
oscillation amplitude of high-energy v^ ^ v., in the presence 
of matter (Liu et al., 1998). While the combination of these 
two analyses disfavors V^ ^ Vj compared with v^ ^ Vx, the 
observation of appearance of the newly created Vx would be a 
definitive answer. Presented here are an update to the NC and 
matter effect analyses and recent attempts to resolve such Vx 
appearance. 



2 Data Analysis 

Super-K is a 50 kt water Cherenkov detector employing 1 1 , 146 
photomultiplier tubes (PMTs) to monitor an internal detec- 
tor (ID) fiducial volume of 22.5 kilotons. Entering and exit- 
ing charged particles are identified by 1885 PMTs in an op- 
tically isolated outer volume (OD). Details of the detector, 
calibrations, and data reduction can be found in (Y. Fukuda 
et al., 1998, 1999). Fufly-contained (FC) events deposit aU 
of their Cherenkov light in the ID while partially-contained 
(PC) events have exiting tracks which deposit some light in 
the OD. The vertex position, number of Cherenkov rings, 
ring directions, and momenta are reconstructed and the par- 
ticle types are identified as "e-like" or "/j-like" for each Che- 
renkov ring. In the current 79.4 kiloton-year (1289 day) FC 
sample, there are 3490 single-ring e-like events, 3346 single- 
ring /7-like events and 3477 multi-ring events. 

Based on 1268 live days, this detector has also collected 



2 



1416 upward through-going muon (UTM) events produced 
by atmospheric neutrino interactions in the surrounding rock. 
We required a minimum track length of 7 m in the inner de- 
tector and an upward muon direction. Downward-going neu- 
trino induced muons cannot be distinguished from the 3 Hz 
of cosmic ray muons. 

2.1 Observing disappearance 

The FC single-ring data sample was analyzed to estabUsh the 
region in (sin^ 29, Am^) oscillation parameter space where v^, 
disappearance is taking place. This data sample is dominated 
by quasi-elastic neutrino interactions, which allow the tag- 
ging of the flavor of the parent neutrino via the flavor of 
the observed outgoing lepton. Data was binned by energy, 
outgoing lepton flavor, and zenith angle of the lepton. The 
distance the neutrino has traveled can be inferred from the 
zenith angle, with those neutrinos coming from above having 
traveled only tens of km, and those coming from below the 
whole diameter of the earth. To provide an expected signal, 
70 hve-years equivalent of Monte-Carlo (MC) data was used. 
As pubUshed by Y. Fukuda et al. (1998), the number of of ob- 
served /i-like events in bins with low energies and long base- 
lines is suppressed compared with the non-oscillated MC ex- 
pectations. 

To evaluate the suitability of possible oscillation param- 
eters, a MC expectation was calculated using the oscilla- 
tion probabiUty appropriate for both two-flavor Vx and 

^ Yv oscillations. A '/} comparison of these MC predic- 
tions and the FC single-ring data yielded allowed regions in 
(sin^26,Am^) parameter space for both types of v^, disap- 
pearance oscillations. For the case of ^ v.s, matter effects 
were included, although at the ~ 1 GeV neutrino energies 
in this data sample the matter effects are small and the two 
resulting allowed regions very similar. This analysis results 
in two sets of neutrino oscillation parameters which are con- 
sistent with disappearance oscillations for a Vx oscillation 
partner or a v^. For details of this fit see Y. Fukuda et al. 
(2000). 

2.2 Neutral Current Sample 

A set of cuts were devised to isolate a comparatively NC rich 
data set. These events must be contained, neutrino-induced 
events with multiple rings, the brightest of which must be 
identified as an electron. Furthermore, to improve the angu- 
lar correlation of the observed particles to their parent neu- 
trino,the total visible energy must be greater than 400 MeV. 
This results in a mean angle difference between the parent 
neutrino and the reconstructed event direction of 33°. When 
tested on the MC data, these cuts produce a data sample con- 
taining a 29% fraction of NC events. The FC single ring 
sample discussed above contains only ^ 6% NC events. Fig- 
ure 1(a) shows the zenith angle distribution of these events 
with predictions from the MC. 

Since Vj do not interact via a NC channel, if the are 
changing into given sufficient distance to do so then the 



long-pathlength zenith angle bins coming from below will 
have comparatively fewer NC events and the ratio of up-to- 
down events will be reduced. However, Vx do interact via the 
NC channel. Thus, if the explanation of the mu-like event 
suppression observed in the FC data is that are changing 
to Vx, then the NC event rate coming up from below will 
remain unchanged, as will the up/down ratio of NC events. 

We define "upward" as a cosine of zenith angle less than 
—0.4 and "downward" as greater than +0.4. There are 465 
upward events and 438 downward events in the NC-enriched 
sample. Figure 1(b) shows the Am^ dependence of the ex- 
pected up/down ratio in the case of full mixing (sin^ 26 = 1). 
For Am^ = 3.2 x 10~^ eV^, the data are consistent with v^u <-> 
Vx, while the up/down ratio predicted by <-> V, oscillations 
is 3.4a too low. 

Using this up-to-down ratio helps to cancel some system- 
atic uncertainties, particularly the rather large uncertainties 
in the NC cross-sections. The total systematic uncertainty 
in the up/down ratio of the data and MC is estimated to be 
±2.9% (Y. Fukuda et al, 2000). 

2.3 Partially Contained Sample 

To look for the presence of matter effects upon the oscilla- 
tion amplitude, we turn to a higher energy event sample as 
matter effects are strongest at high energies. The partially 
contained (PC) event sample is created by neutrino interac- 
tions in the fiducial volume of water which produce an out- 
going particle with a high enough energy to escape the in- 
ner detector. PC events in Super-K are estimated to be 97% 
pure v^ CC and result from parent neutrinos with a mean of 
10 GeV. A further requirement of a visible energy of greater 
than 5 GeV raises this mean neutrino energy to 20 GeV. Fig- 
ure 1(c) shows the zenith angle distribution of these events 
with predictions from MC, as before. Again an up/down ra- 
tio is employed to cancel systematic uncertainties using the 
same angular definition as for the multi-ring sample. There 
are 46 "up" events and 90 "down" events. Figure 1(d) shows 
the Am^ dependence of the expected up/down ratio in the 
case of full mixing. Since matter effects suppress oscillations 
to Vs, the up/down ratio would remain near one in the case of 
^ oscillations. However, for Am^ = 3.2 x 10~^ eV^, 
the observed up/down ratio is consistent with v^ <-> Vx oscil- 
lations and 2.9o too low for v^ ^ v,. The total systematic 
uncertainty in the up/down ratio is estimated to be ±4.1%. 

2.4 Upward Through-going Muon Sample 

An even higher energy sample of neutrinos is observed by 
Super-K and produces upward through-going muon events. 
These neutrinos interact in the rock below Super-K and pro- 
duce an outgoing muon with enough energy to reach and 
completely traverse the detector The neutrinos which pro- 
duce such events have a peak energy of 100 GeV, maximizing 
any possible matter effect suppression of oscillation ampU- 
tudes. Due to the down-going presence of cosmic ray muons, 
a "vertical" (cos 9 < —0.4) to "horizontal" (—0.4 < cos 9 < 
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—0) ratio is used to probe this sample. Figure 1(e) shows 
the zenith angle distribution of these events with predictions. 
Figure 1(f) shows the Am^ dependence of the expected verti- 
cal/horizontal ratio in the case of full mixing. At the point of 
Am^ = 3.2 X 10^^ eV^, the data are consistent with Vx 
oscillation, while the predicted suppression of oscillations 
via the matter effect on <-> oscillations is again not seen 
at a significance of 2.9o. We estimated the total systematic 
uncertainty in the horizontal/vertical flux ratio to be ±3.3%. 

2.5 Evaluating <-> versus <-> 

All four data samples discussed (FC single ring, NC, PC, 
and UTM) are independent. No data found in one is present 
in another. Thus, the results of the four can be statistically 
combined to address the question of versus v^. The FC 
data sample sets the allowable range of (sin2 2e,Am2) that fit 
the observed disappearance oscillations for both Vx and 
Vj. Each pair of oscillation parameters in those allowed re- 
gions is used to generated an expected neutrino signal for 
each of the other three samples. This allows a comparison 
between the up/down or vertical/horizontal ratios observed 
in each of the three samples and that expected for both hy- 
potheses. The data samples being independent, the sum of 
the squares of the three deviations between data and expec- 
tations should follow a distribution with three degrees of 
freedom. A (sin^ 29,Am^) pair is disfavored at the 90(99)% 
C.L. if the value of the sum is greater than 6.3(11.3). For the 
case of v^, ^ v,, all parameters allowed at 99% C.L. by the 
FC single-ring case are disfavored by least the 99% C.L. 
when the other three data sets are examined. However, all 
oscillation hypotheses in the ^ Vx FC 90% C.L. allowed 
region remain allowed by the three additional data sets. 

2.6 Evaluating <-> (Vx -|- V,) 

While the above analysis disfavors pure oscillations, 
Vj could still be present as a sub-dominant partner in domi- 
nant Vx oscillations if Vx and v^ form a close mass dou- 
blet (Fogli et al, 2001). To test the extent of such mixing 
which would produce an expectation which remains com- 
patible with the observed Super-K data, a global was 
formed between aU the atmospheric neutrino data present in 
Super-K (FC, PC, UTM, multi-ring NC, multi-ring ^-like, 
and upward- stopping muons) and MC expectations including 
the appropriate level of x appearance and matter effects. Sys- 
tematic errors were treated in a similar fashion to Y. Fukuda 
etal. (1998). 

This global fit has 190 degrees of freedom. A pure Vx os- 
cillation partner produces expectations which have a x^ dif- 
ference from the data of 175.0 at the best fit point. Pure v^ 
as the oscillation partner yields x^ = 204.8 at its best fit, and 
is disfavored by > 5a compared to pure Vx- If the fraction 
of v^ present in the oscillation sample (sin^ ^) is varied be- 
tween these extremes such that — > (cos^Vx -I- sin^v^), the 
maximum fraction of v.s allowed at the 90(99)% C.L. by the 
complete set of Super-K data is sin^ ^ = 25(35)%. Thus, it is 



still possible for v^ to be present in the atmospheric neutrino 
oscillations as a sub-dominant partner. 

3 X Appearance Searches 

In contrast to looking for (and failing to find) hints of v^, one 
could directly detect the presence of the newly created Vx by 
observing the x appearance associated with a Vx interaction. 
This is difficult, because the x itself has such a short lifetime 
that one must look for its decay products and then deduce that 
a X was their parent. Such events are both rather complicated 
and comparatively rare, as only a small component of the at- 
mospheric neutrino spectrum exceeds the 3.4 GeV x produc- 
tion threshold. Around 75 such interactions are expected to 
be in the existing data set. Additionally, neutrino interactions 
which produce a number of pions or other mesons can mimic 
a X decay within the finite resolution of the Super-K detector. 
However, three separate analyses have looked for x appear- 
ance in Super-K, concentrating on hadronic decay modes. 
All three were developed using the existing atmospheric neu- 
trino MC as background with a x signal injected. To evaluate 
their expected sensitivity, the resulting analyses were also ap- 
plied to a MC atmospheric neutrino sample with x included 
at the level of maximal mixing and Am?' = 3.0 x 10"'', the re- 
maining range of allowable oscillation parameters being con- 
sidered here as a systematic error. 

The first analysis makes use of the standard Super-K event 
reconstruction tools, which use likelihood calculations to com- 
pute to what degree the observed light pattern matches the 
expected distribution for a given set of particles. To isolate 
X events, cuts are imposed on the visible energy, number of 
rings, number of decay electrons, the fraction of energy in 
the brightest ring, the distance from the first ring to the decay 
electron, the max muon momentum, the transverse momen- 
tum, and the particle ID of the brightest ring. 

The second analysis makes use of a neural-net algorithm, 
using variables obtained from the likelihood-based tools as 
input: number of rings found; number of seeds used as the 
input to the ring counting algorithm; number of decay elec- 
trons; particle ID of most energetic ring; number of muons; 
and visible energy. The net was trained and tested on x MC 
data and asked to classify data events as x or not. 

The third analysis makes use of the fact that a x produced 
near threshold will decay symmetrically, while a non-x pion- 
production event will have remaining momentum and skew 
the resulting particles in one direction. Thus, energy flow and 
jet-finding calculations are made in addition to initial cuts 
from the standard likelihood analysis. 

After development and testing using the MC data as in- 
puts, the analyses were then applied to the atmospheric neu- 
trino data set. The events were binned into five bins in cos 0, 
since if Vx are being created via oscillations they will be com- 
ing upward (at long pathlengths). A comparison between the 
data, the MC expectations with no x appearance included, 
and with an expectation including x appearance was made. 
The significances of the potential x excess are around 2a 



4 



Analysis 


Ex 


fit#T 


^exp 


c 


Likelihood 


43% 


66±41,,„,+f^ 


2.0 


1.5 


Neural Net 


51% 


92±35s(a(^23 


2.0 


2.2 


Energy Flow 


32% 


79+'^{stat + sys) 


1.9 


1.8 



Table 1. x-appearance analyses, their efficiencies at saving % events, 
the number of x they observed, and the expected and observed sig- 
nificances of the signals. 



for all three analyses. A fit to the amount of T appearance 
which best describes the data yields the number of observed 
X shown in Tab. 1 . These three analyses are highly correlated, 
so one cannot combine their results in any statistically mean- 
ingful way. While all three are consistent with the presence 
of T appearance, none are statistically significant. 



4 Summary and Conclusion 

Three independent data samples are presented that discrimi- 
nate between the oscillations to either Vx or in the region 
of (sin^ 20,Am^) parameter space favored by the majority of 
the Super-K data. While two-flavor <-> Vj v^, Vx oscil- 
lations both fit the low energy charged current data, the pure 
hypothesis does not fit the higher energy and neutral cur- 
rent samples. Pure <-> Vx neutrino oscillations do fit all 
of the data samples presented, although a simultaneous fit to 
all available Super-K atmospheric neutrino data allows up to 
a 25% (90% C.L.) contribution from v^ as a sub-dominant 
oscillation partner. Additionally, searches for the T appear- 
ance associated with v^ Vx oscillations result in signals 
which are consistent with the presence of x appearance in the 
Super-K detector. 
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Fig. 1. (a,c,e) Zenith angle distributions of atmospheric neutrino 
events satisfying cuts described in the text: (a) multi-ring sample, 
(c) partially contained sample, and (e) upward through-going muon 
sample. The black dots indicate the data and statistical errors. The 
solid line indicates the prediction for ^ v^, and the dashed for 
« V„ with (Am^, sin^2e)=(3.2 x 10"^ eV^,!). The two pre- 
dictions are independently normalized to the number of downward- 
going events for (a) and (c) and the number of horizontal events 
for (e). (b,d,f) Expected value of the corresponding test ratio as a 
function of Am^. The solid horizontal lines indicate the measured 
value from the Super-K data with statistical uncertainty indicated 
by dashed lines. Black dots indicate the prediction for <-» Vx, and 
empty squares for <-> V„ in both cases for maximal mixing. 



